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Abstract: 

In the present paper we have directly computed the Berry curvature terms 
relevant for Graphene in the presence of an inhomogeneous lattice distortion. 
We have employed the generalized Foldy Wouthuysen framework, developed 
by some of us [H H3 [5] . We show that a non-constant lattice distortion leads 
to a valley-orbit coupling which is responsible to a valley-Hall effect. This is 
similar to the valley-Hall effect induced by an electric field proposed in [14] and 
is the analogue of the spin-Hall effect in semiconductors [THl E] . Our general 
expressions for Berry curvature, for the special case of homogeneous distortion, 
reduce to the previously obtained results [13]. We also discuss the Berry phase 
in the quantization of cyclotron motion. 
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1 Introduction 



The latest avatar of crystalline carbon - Graphene [T] - has evolved into a meet- 
ing ground of Particle and Condensed Matter physicists. The reason is that un- 
der certain conditions the electronic excitations in graphene behave like massless 
relativistic fermions moving with (energy independent) Fermi velocity Dp rs 10 6 
m/s, that plays the role of the velocity of light [2]. The honeycomb lattice of 
graphene is made up of two triangular lattices (A and B) and the hexagonal 
Brillouin zone has two distinct and degenerate Dirac points or valleys (K + and 
K-) where the conduction and valence bands meet. In the undoped case the 
chemical potential passes through the Dirac points so that the valence band is 
filled whereas the conduction band is empty. In the Tight Binding Approxima- 
tion, the dynamical equations of motion for the planar electrons (with spin not 
taken in to account) for the above four types can be arranged such that one 
has a four-component Dirac equation at hand. This truly remarkable mapping 
has led to the predictions that a number of Particle Physics phenomena, such 
as Klein tunnelling, Andreev reflection among others, can be (and have been) 
observed in graphene [T]. 

However, from the present literature on graphene, it can be generally ob- 
served that explicit computations in the context of graphene are not performed 
in the four component Dirac form. In the present paper we provide a new ap- 
proach where explicit use has been made of the Dirac framework. We exploit 
a generalized Foldy-Wuthouysen (FW) formalism 3j used for four component 
Dirac Hamiltonian (the reason to use the four component Dirac form is thus 
mainly technical), recently developed by some of us [H [6], where one can 
study the Dirac Hamiltonian with very general type of interaction. The result 
comes as a series in powers of Ti. It is very important to stress that in the present 
work, we have also formulated a method of projection that can isolate the single 
valley effects {via two component Pauli spinors) from the Dirac spinor expres- 
sions and hence our method can be used to study single valley and inter-valley 
interaction effects. 

The method relies heavily on the appearance of intrinsic Berry gauge po- 
tentials [3 [8] that induce a non-commutative geometry (Berry curvature) in the 
electron phase space which subsequently affects the particle dynamics [5] (see 
also |10j for a review). In the FW formalism, in order to properly identify the 
single particle operators (see for example [11]). one needs to perform the same 
unitary transformation U on all the operators (position r, momentum p etc.), 
that (partially) diagonalizes the Hamiltonian H , 

UHU+ = E ; r = UrU + = R + ihUVU + = R + TiA R (1) 

Here E denotes the diagonal energy eigenvalue matrix and A R is the Berry 
potential. One can physically motivate the introduction of r by noting that 
in a Hamiltonian framework it yields the correct expression for the velocity. 
Incidentally the configuration space r (as well as the full phase space) turns 
out to be noncommutative. (The earliest example of a noncommutative space 
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(or phase space) was provided in [H].) This is the "algebraic" reason for the 
emergence of Berry phase. 

On the other hand, intuitively the appearance of Berry phase effects is very 
natural in FW formalism. The Berry potential is induced by the electron spin 
coordinates which are treated as "fast" variables as compared to the position 
coordinates which are considered as " slow" degrees of freedom. The same phe- 
nomena was discovered earlier |13j when molecular dynamics was studied in the 
Born-Oppcnhcimcr approximation. 

The advantage of our method is that one can study all these interactions in a 
unified framework. In fact, for a constant distortion and for zero magnetic field 
we recover the results of [T5] . But obviously more interesting is the new results 
that we obtain for non-constant lattice distortion, which is a new valley-orbit 
coupling due to the spatial variation of the mass gap induced by the lattice 
distortion. This coupling is similar to the valley-orbit coupling in the presence 
of an external electric field and in turn contributes to the valley-Hall effect first 
predicted in [14] . 

In this paper we also consider the semiclassical quantization of the cyclotron 
orbit for an electron in a constant magnetic field and for a constant mass. In 
particular we show the important role played by the valley magnetic moment 
in the determination of the Landau levels. 

The paper is organized as follows: In section 2 we discuss the technicalities, 
construct the (non-commutative) covariant position and momentum operators, 
identify the Berry potentials and compute the semi-classical energy spectrum, in 
a generic framework. In Section 3 we apply the scheme to graphene. In Section 
3 we focus on effects that result from the variable lattice distortion. Section 4 
deals with effects due to external electromagnetic fields. Section 5 is devoted to 
discussions. 

2 Covariant phase space and Berry potentials 

We will explain our formalism through the concrete example of graphene. Around 
the K + and K_ points the energy spectrum is linear and the dispersion relation 
can be approximated as E± (K) = Hvf \K — K±\. In the following vp is set to 
unity so that the effective Hamiltonian near both valleys in the presence of a 
distorted lattice, induced by a real bosonic field A (R) and in the presence of 
an external electric potential V (R) can be written as 

H (P, R) = a.P + /3A (R) - eV (R) . (2) 

where e > so that the charge of the electron is — e. In this expression the 
energy is written in terms of the momentum P = KK. in order that Ti appears 
explicitly during the semiclassical diagonalization procedure described below. 
We have also used the following representation for the Dirac matrices 
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where Oi are Pauli matrices. The dynamical operators P and R have compo- 
nents satisfying the canonical rule [Ri , Pj] = ihSij . We keep the 3+1-dimensional 
representation with the understanding that the components in the z-direction 
are null R z = P z — because graphene is a planar 2-dimensional system. 

Since we are interested in obtaining a semi-classical solution, we first di- 
agonalize the "classical" Hamiltonian |T]) (whose arguments are not operators 

but classical commuting variables noted ^P,R^) by a formal FW unitary trans- 
formation Uq and later reinsert the operators back. The Uq in question is a 
conventional FW matrix, 



, E + A R )+l3a.P 

U (P,RJ = w = - (3) 

/ 2E(E + A (r) ) 



such that the diagonal form of Hamiltonian is, 



U H U+ = e (P,R) = /3-^/p 2 + A (r) 2 - eV (r) . (4) 



where E = P 2 + A (tlj . After this " classical diagonalization" the reintro- 

duction of the quantum operators P and R leads to the following semiclassical 
diagonal Hamiltonian Hp [B], 

e (p, r) ~ £o (p, r) + [[e (p, r) , A*] A p > - [e (p, r) , A*] A R <] + 0(h 2 ). 

(5) 

In the above V+ denotes projection on the diagonal subspace and 



e (p,r)=/5 V /p2 + A(r) 2 - e y(r). (6) 

In this expression the canonical operators (R, P) have been replaced by the 
covariant ones, 

r = R + HA R 

p = P + HA P (7) 

where we have also introduced the so called Berry connections 

A R = V+[A R ] ; A P =V+[A P ], 

defined as the projection on the diagonal of .4 R = % [J/oVpC/q 1 , A F = —i [Uq^rUq] 
Here Uq denotes the matrix Uo (P, R) that is the matrix Uo where now param- 
eters have been replaced by operators. 

Notice the reason for projecting on the diagonal is same as in the conven- 
tional FW case [TT] where this ensures that positive energy solutions do not 
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get mixed up with negative energy solutions. In the present Condensed Matter 
scenario, these correspond to conduction band and valence band. 

From Uq we can now deduce the following covariant coordinate and momen- 
tum 

e x P 

r = R-h y -rE 2 . P = P, (8) 

2E(E + A) > 

where e z is the unit vector in the z— direction and XL = I °^ ^ I . Notice 



-a z 

that the momenta remain unchanged. This will be changed in the presence of 
an external magnetic field. 

From the general expression ([5]) we compute the full diagonal quantum 
Hamiltonian, 

e = /3E- ^e z (p x V r A (r)) E z - eV (r) + 0{h 2 ) (9) 



where now E = yp 2 + A (r) 2 . However, all this is in the 3 + 1-dimensional 
framework and to be applicable to graphene we now need to reduce the model 
dimensionally. This we do in the next section. The coordinate and momentum 
operators r and p are covariant operators because the energy operator ([9]) ex- 
pressed in terms of r and p is gauge independent (i.e. independent of Berry 
connection). 



3 Single valley Hamiltonian for graphene 

In case of graphene one generally considers the two valleys to be isolated and 
studies the electron behavior near one valley. This is achieved by keeping the 
components (1,2) of the Hamiltonian for one valley and (3,4) for the second 
valley. Therefore the 2— bands valley energy can be indexed by r = ±1 such 
that : 

h 



e T =a z \E- r^e z (p x V r A (r)) - eV (r) (10) 

where the eigenvalues ±1 of o z corresponds to the conduction and valence band 
energies of the valley r. The position operator also is now given by the projection 
of |S]) on the components (1, 2), 

r = R, r ; —a, (11) 

2E{E + A) V ' 

Now and for the rest of the paper we consider the properties of an electron in 
a conduction band (the positive energy band of (TIT))) ) of a valley, which means 
that the position operator also has to be projected on this subspace leading to 
the expression 

„ he z x P 

r = R r — (12) 

2E(E + A) v ' 
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From this expression of the valley dependent coordinate operator we deduce that 
the Algebra satisfied by the covariant dynamical variables is a non commutative 
one expressed in terms of Berry curvatures 0^ = d^i A^—d^* Aj where £, rj mean 
either r or p, 

[x, y ] = m 2 ei r v = -m 2 ^ (i3) 

•4-2 

[p if r,-] = -ikS^+ih 2 ®?; = -ih5 ij +T^e jk3 P k V Ri A (14) 

[Px.Pv] = to 2 ©S = (15) 

The Berry curvature Q r J = —j^k obtained here is denoted in the rest of the 
paper simply by 0. We mention that the general expression in (fTS"]) reduces to 
that of [14] who considered valley Hall effect for a constant distortion A. We, on 
the other hand, show that for the more general case of non-constant distortion, 
there are additional Berry curvature components of the mixed form 0fj. In the 
absence of a magnetic field the Berry curvature 0^ is zero, interestingly, the 
energy can be rewritten in terms of the Berry curvature as e T = E — eV (r) + 
7LE(0£r+0^). 

At this level, it is interesting to express the conduction band energy in terms 
of the canonical coordinates P and R which reads 



P 2 + A (R) 2 ^ -Px 

v ' 2E(A + E) 



/2A + E 
\ E 



V R A - eV R V 



-eV(R) 



(16) 

For a constant mass term A = est, we observe a valley-orbit coupling term 
2E{a+e) P x e ^R-^ which is of the same origin as that of the fully relativistic 
spin orbit coupling 2 g(m+_E) q '-(^ > x e ^R^0 m the Dirac equation [4] [15]. In 
the literature this coupling is usually written in the non-relativistic limit only 
j^a.(P x eVR^) because it comes out from the Foldy-Wouthuysen transforma- 
tion which is an expansion in 1/m. The semiclassical diagonalization allows one 
to consider the fully relativistic situation. Interestingly for a non constant mass 
term we observe an additional valley-orbit coupling Kt ( 2 e'Ha+e) ) ^ x ^R.^- 
Therefore a variable mass term mimics an external electric field and can thus 
lead to similar effects as with a true electric field as discuss in next section. 
All this shows that the non-commutative position operator encodes through the 
Berry potential the spin or valley-orbit coupling. In particular, as discussed in 
(4] , the dynamics of the canonical operators turn out to be different from that 
of the non-commutative ones. 



4 Effect of variable lattice distortion A(r) 

As already stated, VrA term induces an additional valley-orbit coupling. More 
precisely, the field A (R) is responsible for a valley-orbit coupling term ~ P x 
VrA similar to the valley-orbit coupling induced by the electric field ~Px VrV^ 
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which in turn is of the same nature as the spin-orbit coupling for Dirac electrons 
in the non-relativistic limit (Pauli term). Our results reveal that a defect of 
the lattice distortion can actually be interpreted as an effective electric field, 
which can be responsible of a mechanical force leading to a valley Hall effect 
as the electric field does [T3]. This is closely related to the spin Hall effect in 
semiconductors [16], [17]). This is a new observation. This clearly reminds us 
of the work of Bernevig and Zhang [18] who showed that shear strain in zinc- 
blende semiconductors such as GaAs induces a spin-orbit coupling term which 
mimics the usual spin-orbit in the presence of an electric field (this mechanism 
in [T5] then leads to a quantum spin Hall effect, which is different from the effect 
considered here ). 

Interestingly within our approach we can study the effect of the valley-orbit 
coupling in both " ultra-relativistic" (A << P) and "non-relativistic" limits 
(A»P). 

In the regime of a large gap A >> P, the energy becomes 
P 2 Tit 

e T « A (R) + - x [3V R A(R) - eV R V(R)} ■ e z (17) 

which shows the same behavior than the usual non-relativistic limit energy of the 
Dirac Hamiltonian. For a constant A, the valley-orbit contribution is cancelled 
by the presence of the mass-gap term 4A 2 which is large compared to electron 
momentum (this is similar to the mass term 4to 2 c 2 in the Pauli term of non- 
relativistic Dirac electron). But this coupling can be greatly enhanced in the 
presence of a very sharp variation of the lattice distortion A. This effect could 
be much larger than the valley-orbit coupling in the presence of an electric field 
if VrA >> eVR,V. Such very localized lattice distortion defect might be built 
by adsorbing the graphene on two different substrates in contact creating in this 
way a line defect. 

In the opposite regime A << P we have 

e T = P- tt^P x [V R A - eWnV] ■ e z (18) 

Here the valley-orbit coupling is independent of the magnitude of the distortion 
field and can also be greatly enhanced by a very sharp variation of the distortion 
A. 

Therefore in both cases, a large gradient of the lattice distortion will generate 
a bigger valley-orbit interaction than with the presence of an electric field. This 
coupling may in principle induce a valley-Hall effect very similar to the valley- 
Hall effect considered in [14], the latter having an electric field and a uniform 
A . 

This valley-Hall effect is better seen in the equations of motion which are 
given in terms of the covariant operators. The above noncommutative algebra 
(|13[) leads to a non-trivial modification in the dynamics of an electron in a 
conduction band of a valley which can be written 

r = V p e{l-hQ pr )-hpx 0, p = -V r e (1 - hQ pr ) (19) 
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where we have defined the curvature vector = Qe z and V p e (1 — TiQ pr ) has 
the following meaning V p j£ (<% — TiOfJ). 

In more details at the first order in h we obtain for the momentum dynamics 



p = -^V r A + eV r ^ (20) 

Ej 



such that the velocity reads: 



t-!-£(Y*4^|V,A-W|x.. (21) 



E 2E 3 VV E 

Importantly, the lattice distortion A (R) induces an anomalous velocity p x 
even in the absence of an electric field. Suppose that the distortion is in the y— 
direction. Then we get a spontaneous (non-steady) current both in the longi- 
tudinal (y— direction) and in the transverse a;— direction due to the anomalous 
velocity. Then electrons will accumulate irrespective of their magnetic moment 
at the longitudinal edge. And a electric potential between the edges in the y— 
direction can be measured. In the same manner, we will have a spontaneous 
separation of the electrons depending on the value of their magnetic moments 
and thus an accumulation in the transverse edges. A steady currents requires 
an electric field. If the mass- valley-orbit coupling is stronger than the electric- 
valley-orbit coupling, then the dominant transverse current is due to the lattice 
distortion and leads a transverse velocity v x ~ ( p + £ A ) V y A which de- 
pends of the valley index but is charge independent. Therefore there will be a 
transverse current in the a;— direction j x (r) ~ e4p J dp g(p, r) ^ V y A 

where g (p, r) is the statistical distribution of carrier, leading to a Valley Hall 
effect (separation of the carrier depending of the valley). Here we considered 
the electron in a ballistic regime neglecting scattering with the impurities. If 
the the system is not connected to the leads along the a;— direction, carriers of 
different valleys accumulate near the edges of the sample as discussed in [14] . If 
the density of up and down magnetic moments are equal there will be not po- 
tential difference in the a;— direction. To observe a transverse electric potential 
one needs to control the chemical potential of both species as suggested in [14] . 

We now consider the case of the graphene in an external electromagnetic 
field. 



5 Graphene in a magnetic field 

The magnetic field interaction is introduced via minimal coupling and the Dirac 
Hamiltonian becomes, 

H (P, R) = a ll + f3A (R) + V (R) (22) 

with II = P+eA the covariant momentum. The semiclassical evolution of a 
Dirac electron in an electromagnetic field for a constant mass term was first 
derived in reference [9]. For the generalization to a variable mass we use the 
general method of Hamiltonian diagonalization of [B] . 
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5.1 Hamiltonian diagonalization 



As shown in [5J, the introduction of the gauge potential is a straightforward 
generalization of Eq. (fTU)) and one obtains for the valley r the following diagonal 
representation of the valley Hamiltonian for a carrier in the conduction band 



r 2 + A (r) 2 - r-^e z (tt x VA (r)) + ehEQB, - eV (r) (23) 

with r = H+hA R , tt = II— ehA R x B and B = V x A. The covariant dynamical 
operators are explicitly given by 

„ he z x II , „ eTill , , 

r = R t ^ — and tt = n + r ; — B z (24) 

2E(E + A) 2E(E + A) v ' 



where E = y 7r 2 + A (r) . We note that in the presence of a magnetic field 
the momentum get a Berry connection a well. Therefore the new commutation 
relations are, 

[r Xl r y ] = m 2 Q = -m 2 ^ (25) 
[m,^] = -iM ij -ieh 2 e il e l i 3 B z -in 2 ^e jk3 IL k \7 Ri A (26) 

and we obtain in addition the non trivial commutation relation between the 
covariant momenta, 

[tt X: TT y ] = -iefiB z + ie 2 h 2 QB 2 z . (27) 

We notice also that only the component of the magnetic field in the z-direction 
(perpendicular to the graphene plane) appears in the energy and in the definition 
of the covariant momentum. This peculiarities is due to the 2-dimensional 
nature of the graphene. 

From the above expressions it is also natural to introduce the magnetization 
vector (similar to the spin) 

m{- K ) = -e%EQ= e ^. (28) 

so that interaction with the magnetic field reads —mB z with the magnetic mo- 
ment m being an intrinsic one, associated to a valley and originated from Berry 
curvature. The magnetic moment which interacts with the magnetic field in the 
z— direction only is an important quantity which allows to make a distinction 
between the carriers of the two valleys as proposed by Xiao et.al. [2]. Here we 
will show that the magnetization is an essential quantity for the computation 
of the energy levels (Landau levels). 
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5.2 Semiclassical quantization of cyclotron orbit 

Thus to simplify the discussion consider a zero electric field V (r) = 0, a constant 
lattice distortion A (R) = est = A and a constant magnetic field B —Be z . 

We start with the general semiclassical equations of motion for electrons in 
magnetic Bloch bands 

de 

f = — hi: x 0, and 7r = — er x B (29) 

OTT 

where e is the energy of a magnetic Bloch band. We have an anomalous velocity 
contribution it x which is the "dual" of the relativistic Lorentz force rxB 
(for a non-uniform B one must add to the second equation a kind of Stern and 
Gerlach term —rriS/B which couples the magnetization to the gradient of the 
magnetic field). 

Combining both equations we have 

1 ' dCA and^^^fexB) (30) 



l-ehB.®\dTt) {1- eKB.®) \dir 

These are the usual semiclassical equations of motion of an electron in a mag- 
netic Bloch band except for the presence of the term 1 — eS.B.0 which as 
B =Be z and = Qe z becomes 1 — ehBQ. For the graphene case we have 
e = e T = yjir 2 + A 2 — mB. Therefore, the equations of motion to the first order 
in H becomes 

l-ehBQ \ . ( l-ehBQ \ 

which are the same than for a relativistic particle in a magnetic field except for 
the additional multiplicative factor 1 — ehBQ. Note that we have 7r^| = so 
that E is a constant of motion, and also 4rr = and 4?? = 0. 

' at at 



We now focus on the non-relativistic case A 2 >> 7r 2 so that e r 
mB and in terms of the coordinates Eq. (|31|) reads 



and 



2A 2 



—5—). » = M— 5— ) (32) 



(l-eB%Q\ (l-eB%Q\ . . 

TT x = -en v B\ I, ir y = eit x B\ I (33) 

Choosing the gauge A x = and A y = BX = B(x — Tjft) we have Tt y — 

p y + eB(x — T-^ziTy) + T^^fB — p y + eBx (note the importance of the Berry 

connection here) , which gives p y — ir y — eBx = thus p y is a constant of motion. 
Or as (1 — ehBQ) / A is a constant of motion we can write 

x = uj 2 (xq — x) (34) 
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which is the equation of a harmonic oscillator with the orbit center xo = 
which is the usual expression except that it is expressed in terms of the covariant 
coordinate and momentum and not in terms of the canonical ones. The cyclotron 
frequency uj — (1 — ehBQ) is also the usual result for the quantization of a 
non relativistic electron in a magnetic field except that it is slightly corrected. 
But this correction is actually negligible for the computation of the energy levels 
and we have 

E n *±(^(n + \- T ^ + ^ + 0(tf) (35) 

where the contribution -| comes from the magnetization term. At the semiclas- 
sical order this expression can also be written 



E n ra±WA 2 + 27LBe^n-|-i-~j + O (h 2 ) (36) 

Therefore the ground state is not degenerate as there is only one possibility to 
realize it n — and magnetization up. We mention here the lattice distortion in 
the present case that leads to a mass in the Dirac equation is imposed by hand. 
In a very interesting paper |19j . an explicit mechanism relying on an electron- 
lattice interaction in the presence of an external magnetic field is proposed for 
the lattice distortion and which leads to the energy levels Eq. (|31)j) . 

It is interesting to show that one can easily recover this result from the 
semiclassical Bohr-Sommcrfcld quantification rule 

P x dX =27rft(n + l/2) (37) 

From J24}, tt = P+eA + r 2E e ™^ A) B, we deduce that P x = rr x - t 
(for A x = 0, and IIj can safely replace by tt x at this order in Ti). In the same 
manner gives dX — dx — r 2 e\e+A) ( as ^ * s a cons t an t OI motion). This 
leads to 

ekB 



{^ ^(E + A) ) = 2heB{n + 1/2) > (38) 

where tt 2 = ttI+tt 2 . In the case A 2 >> 7r 2 , we can replace ir 2 by (E n + mB) 2 - 
A 2 and neglect the contribution t e ^ A j , so that 



E n w ^A 2 + 2hBe (n+^- (39) 

and we retrieve the previous result. Note a similar derivation based on the 
Mikitik and Sharlai's theorem [5D] about the Berry phase in the Onsager semi- 
classical quantization of cyclotron orbits is provided in [21j . 

In the opposite regime A 2 << tt 2 , tt 2 w e 2 and (|38|) therefore becomes 
E 2 — rehB — 2TieB (n + 1/2) so that the Landau levels are given by 

(40) 
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Interestingly from this semiclassical analysis we retrieve the exact result. It is 
important to stress that here the contribution ^ does not come from the mag- 
netization term which is zero here but from the Berry connection in momentum 
j4 p which in the non-relativistic case turns out to be negligible. 

6 Conclusion 

The two-dimensional carbon crystalline honeycomb structure - graphcnc - has 
created an enormous amount of activity because of its striking electronic prop- 
erties. Near each of the two distinct and well separated valleys (or Dirac points) 
in the Brillouin zone, the energy dispersion of electrons has a linear structure 
with electrons moving with an (energy independent) Fermi velocity. It means 
that the electrons can be regarded as massless relativistic two dimensional Dirac 
particles. Thus a number of exotic relativistic phenomena are, in principle, ob- 
servable in graphene. It is possible to simultaneously consider the two valleys 
in a four dimensional framework and this allows us to apply the machinery of 
four component Dirac electron in the planar system of graphene. 

In the present paper we have exploited a generalized form [6j of the well 
known Foldy Wouthuysen [3] transformation to generate the energy spectrum 
(in a semi-classical quantization scheme). We construct the covariant position 
and momentum operators which, on account of the intrinsic Berry potentials, 
satisfy a non-commutative phase space algebra. The induced Berry curvature 
terms directly affect the electron dynamics, as for example by generating an 
anomalous velocity term. 

The main advantages of employing the generalized Foldy Wouthuysen for- 
malism [6] is twofold: 

(i) Different forms of interactions can be studied in a unified way. 

(ii) The approximation scheme leads in a straightforward way to semi-classical 
results which are thus valid for a small ratio of the electron wavelength to the 
typical spatial scale of external perturbations (either in mass or in external po- 
tential) . 

Referring to (i) above, in the present paper we have discussed effects of an ar- 
bitrary position dependent lattice distortion which gives rise to a novel "valley- 
orbit" coupling (similar to the spin-orbit coupling) where one observes that the 
lattice distortion induces an effective electric field. Obviously our general ex- 
pressions (such as Berry curvature) reduce to the well known results of |14j for 
constant lattice defect. It is important to stress that in both cases the intrinsic 
inhomogeneity of the lattice drives the effect without any external electric field. 
We have also studied of effects external electric and magnetic fields. 
In the context of (ii) above, we have shown that one can study both the "ultra- 
relativistic" and "non-relativistic" limits in a single setup depending on the 
effective mass-gap (or lattice defect A) being smaller or larger than the electron 
momentum. 

Acknowledgement 1 The authors acknowledge fruitful discussions with F. 
Piechon and J.N. Fuchs. We thank the referees for the constructive comments. 
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